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Colored tensor models generalize matrix models in arbitrary dimensions yielding a statis- 
tical theory of random higher dimensional topological spaces. They admit a 1/N expansion 
dominated by graphs of spherical topology. The simplest tensor model one can consider 
maps onto a rectangular matrix model with skewed scalings. We analyze this simplest toy 
model and show that it exhibits a family of multi critical points and a novel double scaling 
limit. We show in D = 3 dimensions that only graphs representing spheres contribute in the 
' double scaling limit, and argue that similar results hold for any dimension. 
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I. INTRODUCTION 

Random matrix models [l[ provide a statistical theory of random discretized Riemann surfaces. 
The amplitudes of the ribbon Feynman graphs of their perturbative expansion support a 1/N 
__h expansion 0, 8 (where N is the size of the matrices) indexed by the genus of the surfaces. In the 

large N limit the planar graphs corresponding to surfaces of spherical topology dominate 0, 0] . 
Statistical models of fluctuating geometry can be thought as providing either a regularization or 
a fundamental description of quantum gravity [f|. The large- N limit of matrix models yields 
\Q \ in two dimensions an analytic description of dynamical triangulations 0, Sji whose link to non- 

critical string theory in the continuum limit is well- understood 0]. Higher-dimensional models of 
dynamical triangulations have not been equally successful in providing a sensible continuum limit 



- 1—1 

X 



for quantum gravity 
dynamical triangulations 



although a non-local modification, which goes under the name of causal 



geometry at large scale 



11 [has produced substantial evidence for the emergence of an extended 



The family of planar graphs [3H5J dominating the large N limit of matrix models is summable 
with finite radius of convergence. When the coupling constant approaches a critical value g c , the 
free energy is dominated by graphs with a large number of vertices and exhibits a critical behavior. 



It is in this regime that the system reaches its continuum limit and the critical exponents can be 
evaluated. 

More general matrix models exhibit a complex critical behavior, specific choices of the potential 



leading to multi critical points [151 ]. They correspond to (m + l,m) conformal matter coupled 
to Liouville gravity |15ll. and map on conformal field theories on fixed geometries via the KPZ 
correspondence [16l4l9l|. The contributions to the free energy (that is the partition function for 
connected surfaces) at fixed genus exhibit a critical behavior when the coupling g goes to the same 
fixed g c . The corresponding critical exponents are such that the double scaling limit N — > oo, 
g — > g c , with N(g — g c ) 1+ *™ = k" 1 fixed combines all genera 2CH22] and leads to a well defined 
continuum theory with fixed Newton's constant 0|. 

Random matrices generalize in higher dimensions to random tensors |23l425l|. whose perturbative 



expansion performs a sum over random hi ghe r dimensional geometries. Although some power 



counting estimates have been obtained 2614321] . and the symmetries of tensor models could be 
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analyzed 3^-35] the key to analytical rather than numerical results, namely the 1/N expansion, 
was missing until recently. That situation has changed with the discovery of such a 1/N expansion 



38l | for colored 3914421] random tensors. The amplitude of their graphs supports a 1/N expansion 
indexed by the degree, a positive integer, which plays in higher dimensions the role the genus played 
in two dimensions. The leading order graphs, baptized melonic [43j], triangulate the D-dimensional 
sphere and form a summable series. When the coupling constant approaches its critical value, 
the free energy exhibits a critical behavior and, like in matrix models, the colored tensor models 
reach their continuum limit dominated by triangulations with an infinite number of simplices. The 
entropy exponent of the melonic series, analogous to the string susceptibility 7 s tring = — 1/2 of 
the 1-matrix model for the pure gravity universal class, is 7 me i on s = 1/2, |43j]. This discovery led 
to the possibility of new analytical investigations of dynamical triangulations models and their 
continuum limit in D > 3 dimensions 43j, |44| . Colored random tensors are a promising tractable 
discretization of quantum gravity in three and more dimensions and the subject is developing fast 
45-4(J. The understanding of the leading (melonic) order of colored tensor models allows the 



study of the coupling of statistical systems to random geometries in arbitrary dimension. Thus 
one can prove that, unlike in two dimensions 501. the Ising model on a random lattice in higher 
dimensions does not exhibit a phase transition 51(, while the dually weighted models (introduced 
in [52jl in two dimensions) counting dynamical triangulations with a non trivial measure factor, do 



54|. 



[531 ], A detailed introduction and review of colored tensor models is 

The Schwinger Dyson equations (SDE) of matrix models translate into constrains (satisfied by 
the partition function) which obey the Virasoro algebra 5514571] . A similar line of inquiry can 
be pursued in higher dimensions by integrating all colors but one in a colored tensor model and 
attributing an independent coupling constant to each effective vertex. The algebra of constrains 
is, unsurprisingly, much more involved, but still manageable at leading order 4J|. The detailed 
study of this algebra will lead to a precise understanding of the various critical behaviors and a 
classification of the continuum limits of colored tensor models. 

In its full generality this classification is for now out of reach and one must must content with 
the study of some simplified toy models where a more detailed analysis can be pursued. This 
is done in the present paper by restricting to a sub sector of the tensor models consisting in a 
matrix model with skewed scalings. A tensor T n i n o with D indices can be seen as a N x N ^ 1 
matrix T n i^, with n = n 2 . . . n D . General tensor interactions correspond to arbitrary contractions 
of indices, hence they do not respect this splitting. In the sequel we restrict to the subclass on 
interactions which do, that is we only take into account interactions of the form Tr[(TT^) p ]. In 
order to obtain a sensible theory one needs to adapt appropriately the scalings of the terms with 
the large parameter N. This toy model is natural once one takes into account the full SDE's 
of colored tensor model: it consist in restricting to a model whose leading order SDE's close a 
Virasoro algebra. The matrix model with skewed scalings can also be interpreted as a model of 
random surfaces with patches. 

Using techniques consecrated for matrix models, but taking into account the effect of the unusual 
scaling we present in this paper the full solution of the toy model via the all orders SDE's. We 
generalize for all D the multi critical points of matrix models [TH], and find that their susceptibility 

exponents are 1 , identical with those of multi critical polymers 58]. We the particularize for 

convenience to D = 3 and show that the model admits a double scaling limit N — > 00 g — > g c with 



N(g - g c ) 



fixed. We then go on to show that the double scaling limit we uncover is 



not a summation over topologies: indeed only planar graphs with trivial topology contribute to the 
double scaling limit. Viewed as graphs of a colored tensor model in D = 3, these graphs always 
represent spheres. We argue that similar results hold in all dimensions. 

This paper is organized as follows. In section [TT] we present the colored tensor model and its 
relation to the toy model of a matrix with skewed scalings. In section Hill we analyze the toy model 
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and its multi critical points. In section ITVl we deduce the all orders SDEs and present their iterative 
solution. In section |V] we present the double scaling limit of our toy model in D = 3, and in section 
IVII we discuss the implications of these results for general tensor models. 



II. FROM COLORED TENSOR MODELS TO A MATRIX MODEL WITH SKEWED 

SCALINGS 



Our starting point is the independent identically distributed colored tensor model with one cou- 
pling. We denote n,, for i = 0, . . . , D, the -D-tuple of integers Hi = (nu-i, . . . , rijo, . . . , na+i ), 



with riik = 1, . . . , N. This N is the size of the tensors and the large N limit defined in 36H38I] rep- 
resents the limit of infinite size tensors. We set riij = n^. Let ■i/'t., ipt , with i = 0, . . . , D, be D + 1 
couples of complex conjugated tensors with D indices. The independent identically distributed 



(i.i.d.) colored tensor model in dimension D |54l ] is defined by the partition function 

e ~N°F N (X,X) = Zn(k ~ x) = j ^ # e -S(^) ) 

w vo=ee « + NDi L m e n < + NDi D-w d n % ■ m 

i=0 n n i=0 n i=0 

denotes the sum over all indices riij from 1 to N. The tensor indices nj,- need not be simple 
integers (they can for instance index the Fourier modes of an arbitrary compact Lie group, or 
even of a finite group of large order [5^j ) . The partition function of equation (|2.ip is evaluated by 
colored Feynman graphs 39- 4lt] . The tensors have no symmetry properties under permutations of 



their indices (i.e. all iplt,^. are independent). The colors i of the fields ib % ,ih l induce important 
restrictions on the combinatorics of the graphs. They have two types of vertices, say one of positive 
(involving tp) and one of negative (involving tjj). The lines always join a ip l to a ■0 1 and possess 
a color index, i. Any Feynman graph Q of this model is an orientable simplicial pseudo manifold 



[40l . |42j] and the colored tensor models provide a statistical theory of random triangulations in D 
dimensions, generalizing random matrix models. 

One can picture the topological space associated to a graph in a rather simple manner. The 
vertices of the graph correspond to the D simplices of the simplicial complex. The halflines of a 
vertex represent the D — l simplices bounding a D simplex and have a color. Any lower dimensional 
subsimplex is colored by the colors of the D — l simplices sharing it. In figure [U we sketched the 
dual complex in D = 3 dimensions. The vertices are dual to tetrahedra. A triangle (say 3) is dual 
to a line (of color 3) and separates two tetrahedra. An edge (say common to the triangles 2 and 
3) is dual to a face (and indexed by two colors 2 and 3). A vertex (say common to the triangles 0, 
2 and 3) is indexed by the three colors 0, 2 and 3. 




Figure 1. The dual complex in D = 3. 



The ra-bubbles of the graph are the maximally connected subgraphs made of lines with n fixed 
colors. They are associated to the D — n simplices of the pseudo manifold. For instance, the 



4 



-D-bubbles are the maximally connected subgraphs containing all but one of the colors. They are 
associated to the simplices (vertices) of the pseudo- manifold. 

The tensor indices rijk are preserved along the faces of the graph (which are readily identified 
as bi colored connected subgraphs). The amplitude of a graph with 2p vertices and T faces is [38| 

A(G) = (AAfiV^^ 11 ^ . (2.2) 

We generalize the colored tensor model with one coupling to a model with an infinity of cou- 
plings. First we integrate all colors save one, and second we "free" the couplings of the operators 
in the effective action for the last color. When integrating all colors save one the partition function 
becomes 



b{5 Sym(£°) 

where the sum over B° runs over all connected vacuum graphs with colors I, . . . D (i.e. over all the 
possible -D-bubbles with colors 1, . . . D) and p vertices. The operators Tr^g [tp° , ip°] in the effective 

action for the last color are tensor network operators. Every vertex of B° is decorated by a tensor 
^no or ^no' anc ^ ^ e t ensor indices are contracted as dictated by the graph B°. We denote v, v the 
positive (resp. negative) vertices of B°, and 1 % V ^ the lines (of color i) connecting the positive vertex 
v with the negative vertex v. The operators write 

Tr B #,/]=£( n ^xff n *«) > ( 2 - 4 ) 

n v,vet3° 1=0 lUeB 

where all indices n are summed. Note that, as all vertices in the bubble belong to an unique line 
of a given color, all the indices of the tensors are paired. The scaling with N of an operator can 
be evaluated and the effective action for the last color writes (dropping the index 0) 



with uj(B) a non negative integer 44j, |54|. Attributing to each operator its coupling constant and 

D-l 

rescaling the field to T = i^N 2 f we obtain the partition function of colored tensor model with 
generic potential 

Z = e -N°F(t B ) = j dfdT e -N°-iS{T,T) ^ 

S(f,T) = ^f li T H + z ^t B N^Td^W-^ Tr B [T,T] . (2.6) 

s 

Although in the end we deal with an unique tensor T, the colors are crucial for the definition of 
the tensor network operators in the effective action. The initial vertex of the tensor model described 
a D simplex. The tensor network operators describe (colored) polytopes in D dimensions obtained 
by gluing simplices along all save one of their boundary D — l simplices around a point (dual to the 
bubble B). This is in strict parallel with matrix models, where higher degree interactions represent 
polygons obtained by gluing triangles around a vertex. Each index of a tensor T inherits a color. 
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When evaluating amplitudes of graphs obtained by integrating the last tensor T, the tensor 
network operators act as effective vertices (for instance each comes equipped with its own coupling 
constant). One can represent a Feynman graph either as graphs with D + 1 colors (with the 
subgraphs with colors 1, . . . , D, representing the effective vertices) or directly as a stranded graph 
for the tensor T. The latter is built of stranded lines, and vertices representing the connectivity 
of the indices of the tensors T in an effective vertex. For example, for tensors with three indices, 
both 

T a p l p 2Tp0 p l p 2T b q l q 2T J ,Q q l q 2T c r l r 2TgP r l r 2 , 

T a o x ib 2 TpO p i b 2T p o p i c 2T d o x i c 2T d o q i q 2T a o q i q 2 , (2.7) 

are allowed (leading order uj{B) = 0) vertices. In the representation as stranded graphs one 
represent the index contractions in the vertex by strands. Furthermore, listing the indices of the 
tensors T (respectively T) turning clockwise (respectively anticlockwise) around a vertex, the two 
vertices above are represented in figure [2 




Figure 2. Examples of vertices. 



The full model of eq. (|2,6|) proved for now too complex to be solved analytically. A first step in 
this direction consists in studying simplified models taking into account only subclasses of tensor 
network operators. A first subclass is readily identified: one can restrict to vertices such that the 
tensors share alternatively one and D — 1 indices (the vertex on the left in figure [2]) . Introducing 
the shorthand notation T n i^, where n = (n 2 , . . . n D ), they write 

T n0n f p o n T p o R , ■ ■ ■ = Tr (TTtf (2.8) 

where T denoted the N x N D_1 matrix with entries T n o^. One can check that all such vertices 
have degree uj{B) = 0. The restricted colored tensor model is then a model of a random N x N D ~ l 
rectangular matrix. Note that the passage from a stranded graph representation of a Feynman 
graph to the initial colored representation is trivial: one needs only to collapse all the stranded 
lines into colored lines of color 0, without collapsing the strands in the vertex. In particular, if 
a stranded graph is planar, the associated colored graph is also a planar graph. In particular in 
D = 3 a planar colore d g raph always represents a three spheral]- This is expected to generalize in 



arbitrary dimensions 6l| 



III. A MATRIX MODEL WITH SKEWED SCALINGS. 

The most general matrix model with skewed scalings for a N x A- " 1 matrix T is defined by 
the partition function 

Z = e- N ° F = [ dTdfe- ND - lTlV ^ , V(z) = J2^ P ■ (3-9) 
J P =i 



1 There are several ways to see this: as it is planar the graph has trivial homotopy hence it is a sphere, as it is planar 
it admits a planar jacket [2^, |36|, [3?| and using [H3, 60] it is a sphere, etc. 



ti 



Rectangular matrix models have already been considered in the literature. However, to the 
best of our knowledge, this is the first time that a model in which the two indices n 1 and ft have 
different scaling in N is considered. This different scalings have profound consequences on the 
power counting of graphs (the reader can already infer this by comparing the D dependent scalings 
in eq. (|3.9|) with those of usual matrix models) The rest of this paper is dedicated to solving this 
model. 

The Feynman graphs generated by (13. 9p are ribbon graphs with two kinds of faces: the "heavy" 
faces carrying an index n and the "light" faces, carrying an index nr. They can therefore be in- 
terpreted as random surfaces decorated by patches (corresponding to the light faces). As we will 
see in the sequel, these light patches completely change the continuum limit of the matrix model. 
Note that scalings in N in eq. (13. 9h are the unique scalings which lead to a well defined continuum 
limit. Most of the analysis we perform below relies on classical matrix models techniques. How- 
ever, at almost every step of the way, the results we obtain are very different from their classical 
counterparts. 

The coupling constants t p can be represented as integrals in the complex plane 

t n = -L / du -^V{u) , nt n = -L / du \v'{u) , (3.10) 

2m J c M n+i 2m Jq u n 

where the contour C is the circle at infinity. 

Observables A set of observables is provided by the loop observables, traces of powers of the 
N x N matrix TT^ . They are computed by deriving w.r.t to the couplings, 

1 ( Tr[ (rrt) 1 ) = L s a 

N (1) Zdt q 8tq ^ ' 

We will denote below the generating function of these observables, also known as the resolvent by 
W(z) and its large N limit by Wq{z) 

' Tr 



Remark that the resolvent behaves for \z\ — > oo as ~ for any value of N. The loop observables are 
computed starting from the resolvent via integrations in the complex plane 



i (Tr[(rrt)g]) 

n <rr~ 



— f du u q W{u) . (3.13) 
2m J c 



The free energy Fn is the partition function for connected surfaces. The area of a surface is 
proportional with the number of vertices, that is 



(3-14) 

i=2 

This area can be computed as a single derivative. Renaming the couplings t\ = -, ti = the 
potential can be written as V(z) = M z + z~2i=2 a i zl ) an< ^ t ne derivative of F/v w.r.t. g computes 



( Tr[(MM i ') i j ; 1 f 
gd g F = -N-^U- j\ = ~y~ / du V(u)W(u) , 
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where we used that fact that V(u) is an entire function. It follows that the area of the connected 
surface is A ~ gdg In F (in fact one needs to subtract a finite non universal piece from the free 
energy prior to taking the logarithm). Note that for any analytic function / 

< Tr /(TTt)> = ^-l ' du f(u)W(u) , (3.16) 

hence solving the model consists in determining the resolvent W(u). 



A. Eigenvalue treatment 

For random matrices one can change variable in the functional integral and pass to N eigenvalue 



integrals. This is well known for square matrices and can be done (see [62| and references therein) 
also for rectangular ones. In this case one reduces the integral over T and to the eigenvalues A 
of the square N x N hermitian matrices TT^ . Using [02I ] we can then rewrite the partition function 

as 



N N 

N n-i_ N N D-ly( x ,j 



e~ NDF » = Z N = ([[dXi) J] (A,-A,) 2 II(W 

J ° i=l l<i<j<N i=l 

' d\ ie - ND - lG ^ , (3.17) 



with 

i<j 

The integral in (|3.17p is evaluated by a saddle point whose equations write 



G(Xi) = £ V(Xi) - Y, ^ \ + ^ (E lnA ^-|E ln l A i " A * l) • ( 3 - 18 ) 



dG Tr// , 1 N / 1 2 1 \ . 

In order to solve at leading order (|3.19p . one needs only to note that the last term on the right 
hand side of eq. (|3.19p is suppressed in the large N limit for D > 3. The saddle point equations 
therefore decouple at leading order and the large N free energy is 

Fqo = V(zq) — lnzo , zq the physical solution of zqV'(zo) = nt n ZQ = 1 . (3.20) 



n=l 



The loop observables and resolvent at leading order are therefore 



1 (Tr[(TTt)g]) dFoo / lsdzo q 

hm ! j—. = — = Zn + [ V (Zo) — — = zk , 

N^ooN (1) dt g ° V K °' z J dt q ° 



1 (. Tr _rprp\ 

W (z) = lim ' . ' = — . (3.21) 

1 ^ N^oo N (l) Z - ZQ ^ ' 

In particular Wo(z) has a pole singularity at zq, not a cut singularity. This is the first major 
difference with the usual matrix models, and a direct consequence of the factorization of the saddle 
point equations at leading order. 
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Note that zq has a straightforward interpretation as the leading order connected two point 
function of the model, and the loop observables factor at leading order into two point functions 

1 (Tr[(TTt)<n\ / l (Tr\TT^])\q 

lim 4-± 1 \ . ' V = ( lim i- K . V ) . (3.22) 

The large N factorization of the loop observables can readily be understood in terms of graphs. 

( Tr[(TTt)8]) 

The graphs contributing to the connected correlation ^ are connected vacuum graph 

with a marked vertex (corresponding to the insertion Tr[(TT^) q ]). Chose a tree in any graph 
contributing to this correlation and set its root to be the marked vertex. Going around the tree 
one finds that there exists exactly one contraction which closes a maximal number of faces with 
"heavy" index n. A close inspection reveals that this contraction necessarily connects any two 
consecutive TT^ on the marked vertex via some connected two point graph, and the factorization 
follows. 

Moreover, as zq is the two point function, the eq. zqV'(zq) = 1 is in fact the Schwinger Dyson 
Equation (SDE) 

\ J dTdT^-{T nir , e-^" 1 W(TT)) = Q (3 23) 

combined with the factorization property of the loop observables. Also, zqV'(zo) = 1 is a generating 
equation for trees 

z o = ^ - ^2p f f4 =9-J2 pa P z o ' ( 3 - 24 ) 

11 p=2 11 p=2 

and arbitrary weighted trees exhibit multi critical behaviors. The derivative of free energy F^, at 
leading order can be computed either directly from (|3,20p or from f|3. 15|) and 

QdgF^ = I du V(u)—^— = -V(z ) . (3.25) 

2m J c u- z 



B. Multi Critical Points 



Inspired by [15J we choose for potential a polynomial of degree m. The saddle point equations 



zoV'(z ) = 1 5 = z + ^2 Q a Q z o > ( 3 - 26 ) 

prove that the coupling g is a polynomial in zq. The tensor model achieves its continuum limit 
when 



= , (3.27) 

OZ 



and a m'th multi critical point when 



da d m ~ 1 a d m a 

dzo (dzo)™- 1 (dz ) m T y ' 

A minimal realization for a m'th multi critical model is then obtained if the saddle point equation 
is 

g = g c - (z c - z ) m , 9c = z? , (3.29) 
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corresponding to a choice of potential 

m— 1 ., -, m— 1 



lJ 9 c ±* g+i 



m— q— 1 



9+1 



771 



(3.30) 



g=0 * 3 q=0 

where we recall that V(0) = and the coefficient of z is -. The leading order resolvent is then 



1 



(Zc ~ ZO) 



E {g c - g) m 



Z — Zn * — ' (z — Z c ) k+1 ' (z — Z r ) 

u k=0 y CJ k=0 y cl 



(3.31) 



and the leading non analytic behavior of the derivative of the free energy when g — > g c computes 



gOgFoo 



{g c -g) m F OQ = (g c - g) 1+ ™ 



(3.32) 



corresponding to a susceptibility exponent 7 m = 1 — 1/m, identical with the one of the multi 
critical polymers of (58|. This is not surprising, taking into account the earlier remark on trees. 
The g — > g c limit is a continuum limit as the average area, A ~ d g InF ~ j—^ ; , diverges. 



IV. SCHWINGER DYSON EQUATIONS AT ALL ORDERS 

In this section we go beyond the leading order and derive the SDE's at all orders. We sub- 
sequently present the iterative solution order by order in N. In order to fit the scalings in N of 
various terms one needs to take specific choices for the dimension D. We will do this in the next 
section. 

The derivation of SDE's follows the classical path of 55|-57|, up to the unusual scalings of 
various terms with N. We reproduce it below for completeness. For q > 1 we write a SDE 



dT- 



5T n 



q-l 



r=l 



N( Tr[(T^T) q ]\ + ( Tr[(TTt)«- r ] Tr[(T^T) r 



N 



D-l 



Tr[(TT t ) ? 



-N D - 1 Y,Pt P { Tr[(TT^) q+p ]\ = , 



where the double trace term is absent if q = 1. The case q = is special and leads to 



N D (1 



N 



D-l 



-£pt p (Tr[(TT^}) = 



(4.33) 



(4.34) 



The SDEs can be written as 
L m Z = , 



L 



-D+2 



+ 1) 



-(N 
Li = -(N 

L = N D + Y J Ph 

p=i 



d_ 

dta 



-2D+2 







dt 

r=l * 



+ Y1 &p 







p=i 



dt 



Q+P 



dt! 
d 



P =i 



'dt 



p+i 



Of 



q+p 



(4.35) 
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and, like for the usual matrix models, the L m 's respect the Virasoro algebra. Introducing the 
derivative of the resolvent 



p=i 



(4.36) 



a straightforward computations expresses the double trace observable as 



1 ^(Tr[(TTt)*-]TV[(rrtr 

2^1 yt+l Z_> 



t=2 r=l 



AT^Wfo z) + iV 2 (W(z) 2 - 2W{z) + -) . 



(4.37) 



Adding the equations (|4.33p and (|4.34p with well chosen coefficients leads to 



1 i 1 {Tt[(TT^) 



! V Tr[(TTt)]) x x ( T r[(TTt)^; 



p=i 



(TV[(rrt)^]) ! / Tr[(TT t ) 9 ~ T '] Tr [(T t T) r 

q=2 { 1 ) r=l 



1 ( T r[(rrt)^; 



p=i 

Recalling that 
we obtain 



. 



(4.38) 



pt p 



2ni 



I 1 1 ( Tr[(TTt)g] 

9=1 



du —V'(u) 

c uP 



(4.39) 



Tr[(rrt)«; 



9=1 



i ( Trf(Trt)9- r l Tr[(TtT) r ] 



g =2 r=l 



4 y J_ / dn ±y (l 



1 ( T¥[(TTt)9+P] 



Z 9+1 



o, 



(4.40) 



and a straightforward computations finally yields the loop insertion equation 



1 



uV'(u) , . 
du y -^W{u) 



2m J(j z — u 

= (l- zV'{z)^j W{z) + N- D+2 (zW{zf - W(zfj + N~ 2D+2 W (z, z) . 



(4.41) 
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Note that for D = 2 we obtain the usual loop insertion equation for a non hermitian matrix (see 
for instance [55J]). The loop insertion equation simplifies further. Indeed, at leading order Wq{z) 
depends on ti only implicitly trough zq. This generalizes order by order to W(z). Noting that 



dz 



Zo V'(z )-l = 0^ zV'(z) dzo + ^nz^dtn = 0^ -± 

Ota 



q4 



n=l 



zV'(z) 



l , (4-42) 



z=z 



where we supposed that zV (z) ^ (which holds for instance for multi critical models), we 

L J z=z 

can evaluate 



q =i z m i q =i z \zV'(z) 



z=z 



and for z > zq we obtain 



W(z,z) 



zz Q d Zo W(z 







i 




(z - z )2 


zV'(z) 




Z = ZQ 



and the equation (14.41 H becomes 



uV'{u) 
2lTl J c z — u 



da 



W(u) = (l - zV'(z)j W(z) + N~ D+2 (zW{z) 2 - W(z] 

zzo d Zo W{z) 



-2D+2 







1 




(z - z ) 2 


zV'{z) 




Z=ZQ 



(4.43) 



(4.44) 



(4.45) 



A. Iterative Solution 



The loop insertion equation (14. 45ft can be solved order by order in N. We start by expanding 
the resolvent in powers of N, 



W(z)=W (z)+J2 N ~ nW n( Z ) ■ 



(4.46) 



From eq. (j4.45H we conclude that the first non trivial correction is Wd-2> and it respects the 
equation 



uV'(u) 
2m J c z — u 



Note that 



da. 



W D - 2 (u) - (l - zV'{z])w D - 2 {z) = zW 2 - W 



zq 



(z - z Q f 



du 



uV'{u) 1 



2m J c z — u (u — zo) n+ 

n 



E 

p=0 
n-l 

E 

p=0 



1 



1 



(z - z ) p+1 {n-p)\ 



1 



1 



(z - z )P+ 1 (n-p)\l 



uV'{u) 
uV'{u) 



(n-p) 



(z - Z ) n+1 
1 



z (Z - Z ) n + l 

(n-p) 

i 

20 



(4.47) 



(4.48) 
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where we took into account ZqV'(zq) = 1. Proposing the ansatz 

a b 



W D . 



+ 



(z - Zq) 3 (Z - Z ) 2 ' 



The equation writes 
1 



.(z-z )H 
with solution 



uV'{u) 



+ 



1 1 



zo 



z — zq2 . 



uV'{u) 



zq/ \ Z — Zq 



uV'{u) 



zq 



zj (z - z )' 



(4.49) 



(4.50) 



zq 




i 


uV'{u) 









ZQ 



uV'(u) 



2{|W'(u)l' ) 2 

J ZQ ) 



(4.51) 



The next non trivial term is Wzd-4- Subsequent terms are W3D-6 an d W2D-2 (generated by 
the last line in eq. (|4.45p ). We obtain non trivial corrections W n for n = p(D — 2) + q(2D — 2). In 
order to derive the full solution we will in the sequel chose D = 3. 



V. THE DOUBLE SCALING LIMIT IN D = 3 

In D = 3 we obtain a nonzero correction W n for all n 

W{Z) =W (Z) +J2 N ~ n Wn{z) ■ 



(5.52) 



i=i 



Note that for any even D (both D = 4 and in D = 2 for instance) nontrivial corrections W n arise 
only for even n. The equation (|4.45p translates for W n into 



J_ f du HXM Wn{u) . (i- z v'(z))w n (z) 

2,111 J c Z — U \ J 

n-l 

= z^ W p (z)W n - 1 - p (z) - W n -i(z) + f w 

p=Q (z — Zq) 2 uV'{u) 



ZZq 



—d Zo W n -±{z) 



(5.53) 



where by convention W„ = for q < 0. To solve the equation for W n , we propose the ansatz 



*=1 

and substituting we obtain for n > 1 

2n p-l 1 1 



z - Z 



(5.54) 



p=l 



(p-fc) 

20 



fc=0 

n— 1 n— 1 

- Zq) 

W p( Z ) W n-l- P (z) + Z Y J W p {z)W n - 1 - p {z) -W n -t( 
p=0 p=0 



zq 



(z - zq) uV'{u) 



—d Zo W n - A (z) + 



-0 



z» 



(z-z ) 2 uV'{u) 



—d Zo W n -4:(z) 



(5.55) 
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Exchanging the sums on the left hand side we obtain 

2n-l 2n 



E 



1 



(z - z ) k + 

k=0 y u/ p=k+l 



(p-fc) 

2(1 



(p - A;)! 

n—1 n—1 
Z ^ W p (z)W n ^ p (z) + {Z- Z ) ^ Wp^JWn-i-pCz) - W n -l(i 
p=0 p=0 



(z - z ) uV'{u) 



-d Z0 W n . 4 (z) + 



(z-z ) 2 uV'(ii) 



— 9 Z0 VF n _ 4 (z) 



(5.56) 



The equation ()5.56|) allows to determine iteratively the coefficients fn(zo). In order to prove 
that a multi critical model admits a double scaling limit we are interested only in the behavior of 
these coefficients when zo approaches z c . Using the solution at first order we see that at leading 
order in z c — zq, 



f?(zo) 



d\ 



(z c - Zo)™" 1 ' 



fh*>) 



(z c - ZqY 



(5.57) 



with d\ and d\ some numerical coefficients. We will show by induction that the generic scaling 
formula 



mz ) 



(z c - z ) n ( m + 1 )- k 



(5.58) 



holds. The scaling at order n is fixed from (|5.56p . We analyze one by one the scalings of the various 
terms on the right hand side for n > 2 
The first term we analyze is 



n-l 



n-2 



^ W p {z)W n ^ p (z) = 2W W n . l + W^Wn-x-piz) , 

p=0 p=0 
where the second sum is absent when n = 2. It writes 

n-2 2p 2n-2-2p fk\( „ \ f k 2 ( y \ 
2n-2 9 -t / \ n-2 2p 2n-2 

= £ +EEE 77^fZ'MA-U*) 



(5.59) 



E 



p=0fei=l t=2 v uy 
1 1 

(Z - Z ) fc+1 (Z C - 2 )(n-l)(m+l)-fc+l ' 



(5.60) 



hence the first line of the right hand side of (|5.56p is dominated when zo — > z c by the first term. 
The derivatives of W n _4 behave for zq —> z c like 



d Z0 W n - 4 ~ ^ 



1 



1 



(Z-Z ) fc+1 (Z C - Zo )("-4)(m+l)- 



•fc+1 



hence the second line on the right hand side of f|5.56j) scales at most like 

1 1 



E 



(z-z ) fc+1 (z c -z )( n - 3 )( m+1 )- 



fc+3 ' 



(5.61) 



(5.62) 
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and when zq — > z c is dominated by the terms in the first line. Identifying the coefficients of ^ z _^k+i 
in eq. (|5.56|) we conclude that 



2/i 



E f P n^)[uV']^ ~ - _ (ra _ 1)(m+1) „ fc+1 , (5.63) 

p=k+l ^ 6 U ' 

and taking into account that [uF']^"^ ~ (z c — zo) m ~ p+fc i we get 

fn( z o) ~ 7 , (5-64) 

reproducing (|5.58j) . It follows that the leading contributions to the resolvent write 



z ) k d k n 



Z ~ Z ^ ^ (-2 " Z 0) k+1 (Zc ~ ^ ) n(m+1) 

= + § ^ - *° ? d " ^ - z °) m+i ] n ' (5,65) 

Taking into account that z c — zo = (<?c ~~ the resolvent writes when g ^ g c as 

W(z) = -L- + £ (z _^ )fc+1 ( gc - (5.66) 

with k _1 = N(g c — go) 1+ ™. The graphs contributing to the double scaling limit can be easily 
understood in terms of stranded graphs starting from the SDE in equation (|4.4ip . The term on the 
left hand side represents the connection of two effective vertices by a tree line. The terms on the left 
hand side represent a marked vertex decorated by a loop line. The loop line divides the marked ver- 
tex into two smaller effective vertices, corresponding to 7^y( Tr[(TT*) q ~ r ] Tr[(T , 'T') r ]). Either these 

two smaller vertices do not reconnect again (corresponding to the terms involving W(z) and W{z) 2 ), 
i.e. the correlation splits into two connected correlations ( Tv[(TT^) q ~ r \j ^ Tv[(T^T) r ]^ c , in which 
case the total graph is planar. Or the two smaller vertices reconnect with at least a line (the terms 
generated by W(z, z)), in which case the correlation does not split < Tv[(TT^) q - r ] Tr[(T*T) r ]) c , and 
the initial loop line and the new line cross yielding a non planar graph. All the contributions of 
W(z,z) are suppressed in the double scaling limit with respect to those of W 2 (z). It ensues that 
only planar graphs corresponding to spheres in D = 3 dimension contribute in the double scaling 
limit. 



VI. CONCLUSION 

We have shown that a N x N®^ 1 matrix model admits multi critical points and a double 
scaling limit. The model constitutes a toy model for higher dimensional tensor models. In D = 3 
the double scaling limit only graphs having the topology of a sphere contribute. This is expected 
to hold for all D. Indeed, form eq. (|4.41j) we see that the term W(z, z) is increasingly suppressed 
with respect to W 2 (z), hence only planar graphs should dominate for all D. Such graphs are 
homotopically trivial in all dimensions, and are expected to represent spheres. 

Although the model we analyze is very simple, and the (prospective) double scaling limit of 
colored tensor models is expected to be much more involved one can draw several conclusions 
about the major features of such an expansion. First, one should not expect all terms in the 1/N 



15 



expansion to contribute: unlike in usual matrix models, even in the double scaling many terms 
are suppressed. In fact, somewhat surprisingly, one should not expect the double scaling limit to 
be a summation over topologies: it seems possible that only spherical graphs contribute. Third, 
the exact interpretation of this double scaling limit is yet unclear. In two dimensions one notices 
that the free energy in the double scaling limit is a series in k x , where x is the Euler character 
of a graph. As the latter is just the evaluation of the Einstein Hilbert action, one concludes that 
the double scaling limit corresponds to a continuum limit with a finite (renormalized) Newton's 
constant. In higher dimension the precise interpretation of the parameter n is far from clear. 
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